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We consider the orientational instabilities, both homogeneous and spatially periodic, developing 
in a nematic liquid crystal under rectilinear oscillatory Couette flow for director alignment per- 
pendicular to the flow plane. Using numerical and analytical approaches we determine the critical 
amplitude of oscillatory flow instabilities, the critical wave number and the symmetry of the desta- 
bilizing mode. It was found, that by varying of the oscillatory flow frequency the instability changes 
its temporal symmetry. The mechanism of this transition is coupled with the inertia of the nematic 
fluid. We also shown that an electric field applied to the nematic layer can induce switching between 
instabilities with different spatial and temporal symmetries. The complete phase diagram of the 
flow instabilities is presented. 
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I ' I. INTRODUCTION 

Nematic liquid crystals (nematics) represent one of the simplest anisotropic fluid [Q. In order to describe the 
O , macroscopic dynamics of these materials made up from elongated molecules one introduces as additional hydrodynamic 
variable a unit vector n, the director, which describes the axis of local preffered molecular orientation. The strong 
coupling between the director and velocity field provides a number of interesting phenomena in nematics under flow of 
J> i different types, such as pattern formation |||| , director tumbling 0] , autowaves || etc.. The flow behavior of a nematic 
0^ ' strongly depends on the geometry and material parameters, in particular the viscosity coefficients cti . . . ol% known as 
Leslie coefficients. Two of these coefficients, ct2 and 0:3, appear in the director equation and play an particular role in 
hydrodynamic instabilities. c%2 is supposed to be negative for materials made of rod-like molecules, whereas the sign 
of ct2 is not restricted. If the nematic is initially oriented within the flow plane (= x, z plane) then, in materials with 
£SJ ' CV3 < under steady flow (plane Couette or Poiseuille type) the director tends to align at an angle Ofi = tan -1 \J a^j 0x2 
to the flow direction in the flow plane (= x direction), which corresponds to stationary homogeneous solution of the 
■ nematodynamic equations 0]. Such nematics are called flow- aligning. The solution Ofi with appropriate sign of the 
square root was shown to be linearly stable in the case of steady flows . In oscillatory flows the director tends 
to oscillate between ±0//. With a nonlinear velocity profile, as in Poiseuille flow, this solution can lose stability §,§. 
In the case «3 > under steady flow one has tumbling with a director profile depending strongly on the shear rate 
p|. Nematics with 03 > are called non- flow- aligning. The picture of orientational transitions becomes even richer 
if the director is initially oriented perpendicular to the flow plane (y direction). This orientation corresponds, for 
, symmetry reasons, to a stationary homogeneous solution of dynamic equations of nematics. With appropriate surface 



alignment it can lose stability at a certain magnitude of the shear rate via a stationary bifurcation |L0 11 1. In the 
absence of external fields in steady Couette as well as Poiseuille flows the director profile then deforms homogeneously 
in space. In steady Couette flow a magnetic field applied along the initial director orientation changes the type of 
instability from homogeneous to spatially periodic at a critical field strength p0|. In Poiseuille flow the magnetic 



field does not affect the type of instability |LlJ , but rolls can be observed above a secondary instabili ty [111. Classical 
experiments of Pieranski and Guyon (PG) on oscillatory Couette flow ]To| and recent experiments |12| , |13|| show the 
appearance of a roll instability with the critical shear rate depending on flow frequency. It was shown by PG, that 
the temporal symmetry of the instability is influenced by an applied electric field (lO) . A theoretical treatment of the 
problem @|j|]|§ gives the mechanisms of the instability and is in a good agreement with experimental results, 
including the change in symmetry. For oscillatory Poiseuille flow the situation is similar |l6| Jl8| . In all the theoretical 
treatments @,H-||jrj the inertia of the fluid is neglected, which is reasonable for the low frequencies used. In the 



absence of an electric field one then has no changes in the stability scenario as a function of the flow frequency. 

In the present work we investigate numerically as well as analytically orientational instabilities, both homogeneous 
and spatially periodic, induced by plane, oscillatory Couette flow in nematics in the geometry with the director 
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pre-aligned perpendicular to the flow plane in a wide range of the flow frequencies. We look for homogeneous and 
spatially periodic instabilities and study the influence of an electric field on flow instabilities for both, negative and 
positive, dielectric anisotropy of nematic. Recent experiments fl^,|l3| show that high-frequency range could also be 
investigated experimentally as well. 

II. GOVERNING EQUATIONS AND SYMMETRY 

Let us consider a nematic layer filling the space — d/2 < z < d/2 bounded by infinite parallel plates, one of which 
(upper) oscillating harmonically along x with amplitude A and circular frequency lu = 2ixf. The director is oriented 
(rigidly anchored) at the boundaries along y, i.e., perpendicular to the flow plane. An electric field of strength E can 
be applied along z. 

Before any instability develops (sufficiently small amplitude A) the director remains undistorted and the nematic 
behaves like an isotropic fluid with viscosity 773 = a^/2. For further consideration let us introduce the dimensionless 
variables 

z — z/d, t — tu, a = A/d. 
In following we omit the tildes. Then the subcritical velocity profile can be found from Navier-Stokes equation 

0. (1) 



~ a cos t , 



z=-l/2 



= 1/2 

Here e v = 1/(ujt v ) with t v — pd 2 /r\ 3 the viscous relaxation time. Thus, our basic state is 

n° = (0,1,0), v° = [«£(*,<), 0,0]. (2) 

Below it will be shown that for the stability problem the inertia of fluid in the basic velocity field can be neglected, 
so that the solution of (Q) is simply v^, — a{z + 1/2) cost. In order to investigate the stability of (||) we linearize the 
nematodynamic equations Q around this state 

n = (0, 1,0) + (n x ,n y ,n z ), 
v = (v°,Q,0) + (v x ,v v ,v z ), 

where the perturbations ni — ni(y, z, t),Vi — Vi(y, z, t), {i — x, y, z} are small. Here we neglect the x dependence of 
perturbations assuming the wave vector of instability if not zero is perpendicular to the shear direction. 

Eliminating v y by use of the incompressibility condition v y = —v z e and the pressure by cross differentiating the 
nematodynamic equations one obtains 



(e v 1 L 1 + L 2 )v x = -e v 1 d yy v z a cos(t) + a' 2 d yyyt n x + 
+ ° 5 2" 2 acos(t)d yyy n z ] 



(3) 



(e" 1 ^ + L±)v z = Fw z + a cos(t)F 2 n x ; (4) 

(G 1 +e d G 2 )n x = (1 - A) _1 [n 2 acos(t) + v x y ]; (5) 

(G 3 + e d Gi)n z = (1 - A) _1 [An XiH acos(<) + (d yy - Xd zz )v z ], (6) 
where Li, Gi and Fi are linear differential operators 



u 


— dyyt , 


L 2 = 


u 


= Oyyt + d ZZ t, 


u = 


Fi 


— Ct 2 dyyyt OL^3 ZZ yt, 


F 2 = 


Gi 


= d u 


G 2 = 


G 3 


= dyt, 


G 4 = 



-{i 2 dl+d ZZ yy), 

+ K + V[ + V'2)d ZZ yy + V[d% 



a R —a, 



03 _ ?k+< d . 



2 w y 2 w v zz > 
-(k 3 d yy + k 2 d zz ), 

-[k 3 dl + k 3 sgn(e a )n 2 E 2 d y + d zzy ] 
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and 



e d = l/{T d u) , r d = 7 irf 2 /^ii , E = E/E F , 

71 = a 3 - a 2 , A = a 3 /o!2 , 
m = ("3 + «4 + a 6 )/2 , r\2 = ("4 + "5 - aa)/2 ■ 



Here fcj = Ka/Kii with if^ the orientation elastic constants, the director relaxation time, £7? = n /d^J ^n/(eol e o|) 
the Freedericksz field, e a the anisotropy of the dielectric permeativity. Notation / j = df jdi has been used throughout. 
We introduce the quantity P p — t v /t c i — pK 11/(77371), which is of order 10 -6 for ordinary nematics. 
The boundary conditions {fully rigid) read 



z=±l/2 
*=±l/2 



0, v z 
■0, n : 



z=±l/2 
z=±l/2 



: 0. 



z=±l/2 



0. 



(7) 



III. STABILITY ANALYSIS 



A. General formulation 



Let us collect all perturbations in the vector Y = (v x ,v z ,n x ,n z ) . From the Floquet theorem one can write the 
model solutions of (|3l-|6|) in the form 



e °t e iqy 



E y^y ku i ( 8 ) 



k— — c 



where a is the growth rate, q is a wave number of instability. Instead of solving the Floquet problem (g) directly we 
will employ a Galerkin expansion E3|: 

00 

Yk = E C nk<f>n(z) , (9) 
n=l 

where c„fc are constant coefficients and 4> n (z) is a set of trial functions, which satisfy boundary conditions @), 
Chebyshev functions ^(J 

, , s _ / T 2m (2z) -T {2z), n=2m-l, . . 

0n[ ' ~ \ T 2m+ i(2z) - Ti(2z), 71 = 2m, liUj 

for u x , n x and n z and Chadrasekhar functions pl| 

, c»\ _ / cosh(A m z)/ cosh(A m /2) - cos(A m z)/ cos(A m /2), n = 2m - 1, , . 

' 1 sinh(z/ m z)/ sinh(f m /2) — sin(^ m z)/ sin(^ m /2), n = 2m ^ ' 

for v z . Here A m and v m are the roots of the appropriate characteristic equations . 

Due to the structure of Eqs.(|[]|) the model solutions have the definite parity in z and t (parity of Fourier 

modes). Thus, the symmetry types are divided into four groups (Table Q). The z-parity of v y is opposite to that 
of v z whereas the i-parity is the same. Type I corresponds to the Y"-symmetry and type II to the Z-symmetry in 
the notations of PG |lQ| . Type I and III suppose that v z can have nonzero time averaged and rolls in this case are 
overturning, i.e. they rotate steadily. Rolls of II and IV type are oscillating. 

Substituting (||)-(||) into (|^-|j) and projecting jl9| onto the Galerkin modes one obtains a linear algebraic system 
for the expansion coefficients c n k of the form 

(A + aB)c = erCc. (12) 
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The expansions (@)-@ are then truncated and the eigenvalue problem ( |l2| ) can be solved. The condition Si[cr(a, q)] = 
for the largest sigma yields the neutral curve ao(q). If one assumes a stationary bifurcation, i.e. 3(cr) = at the 
threshold, then one has to solve the eigenvalue problem 

(A _1 B)c = — -c, (13) 
a 

which is much more convenient. The critical value a c is given by a c = min g ao(q), which also yields the critical wave 
number q c . Using (|l2|) we have checked that the first bifurcation is really always stationary [E2|. 



B. Results without electric field 

The calculations were done for material MBBA, see App. ^ for material parameters. 

For uiTd < 35 the first instability developing in the absence of external fields is a homogeneous one (q c = 0) [10|. 
This could be expected already from the fact that in steady shear flow the homogeneous instability occurs first [ 10 1 . 
At LUTd ~ 35 PG have observed the transition to rolls of type I. Up to now the frequency region lot^ > 35 remains 
experimentally poorly investigated. Our analysis for lot^ > 35 shows that roll instability of type I persists up to 
LUTd ~ 4.1 ■ 10 4 , where a transition to the rolls of II-type takes place (Fig.|l|). At the transition point the wave number 
jumps from 12 to 1.2. We have checked numerically that the appearance of this transition is related to the inertia of 
the nematic described by the parameter e v . By neglecting the inertia term in the Navier-Stokes Eqs. (|^-^|), i.e. for 
e^ 1 — ► 0, no such transition was found. Numerical analysis also shows that thresholds for instabilities of III and IV 
types are higher for all frequencies investigated. 

In order to understand the mechanism leading to this transition let us analyze Eqs.(|J)-(^|). By acting with d y on 
Eq. (^) and with d yy — Xd zz on (Q) one may eliminate the velocities. In the remaining equations for n x and n z we 
use a one-mode spatial approximation 

n x = cos(iTz)cos(qy)f(t), 
n z — cos(nz) cos(qy)g(t), 

appropriate for the symmetry of type I or II. Projecting the equations one obtains coupled second order ODEs for 
f(t) and g(t). Further simplifications can be made by neglecting the terms proportional to P p and e~ 3 . One then 
arrives at 

/ + edmi f + (e^ 1 m^a sini — TO2<2 cost)g = 0, (14) 
g + £<ihg — ha costf = 0. (15) 

Coefficients m,i(q), h(q) are presented in Appendix |A|. In looking for expressions for the critical amplitude we assume 

/ = fosc, g = g s t + g osc for type I, 

/ = fst + lose, g = g osc for type II, 

that corresponds to the separation of the functions / and g into the stationary (index st) and oscillatory (osc) 
parts. Substituting ( |l6|) into (14)-(p"5|) and using time averaging procedure one obtains the critical amplitudes for roll 



instabilities of I and II type 

[m 1 m 2 + (€ v e d )- 1 m 3 ] l 2 



a" = 



[m 2 h-(e v e d ) 1 m 3 ] l 2 



1/2 

1/2' (17) 



Let us analyze the formulas (p"7|). From ( |Al[ ) one has that fa(q) is positive for large q (> 2) and very small q (<c 1). 
Taking into account that h(q) is positive for all q and the expression \m\m2 +rri3(e v ed)~ 1 ] remains positive up to very 
large value of wr^ (~ 10 6 ) one sees from first Eq. (|l7|) that in order to have a threshold 1% has to be positive. Thus, 
only large q (rolls) and homogeneous instabilities can develop with symmetry of I type. In the second Eq.([l7]) the 
expression [m^i — ^(e^ed) -1 ] in the denominator changes sign with increasing frequency from positive to negative 
at the frequency given by 



1 TO 2^1 1 



(ur d )i r = (18) 
m 3 P p 
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A lower limit of the transition frequency is given by setting q = in (|18|). Below (u>Td)t r one has [m 2 h — m^tvtd)^ 1 ] > 
0. Since mi > for all q one needs I2 > to have a threshold, i.e., only the large q instability (rolls) can exist. For 
LUTd > (ujTd)tr one needs negative I2, which provides an instability with small q. One can conclude from ( |l7| ) that only 
the product e v ed plays an essential role in the transition from rolls of I type symmetry to rolls of II type, not e v itself. 

Let us now return to the question of neglecting the inertia term in Eq.(Q) for the basic velocity. Assuming e^ 1 to 
be small (for ioTd <C 1/-P P ~ 10 6 ) one can find the expression for basic velocity v® as an expansion in e~ l . Keeping 
only the first correction changes the threshold amplitude [and, consequently, {ujTd)tr\ by less than 1%, so that e" 1 can 
really be neglected in (|l]). 

The results for the approximate expressions ([l7|) are included in Fig]l|. The formulas ( |l7j ) have an error of about 
30% at small Lord, which grows further with LUTd- This is an obvious consequence of the one-mode approximation. 
Nevertheless, the analytic expressions show all the important features of the numerically calculated thresholds, in- 
cluding the switching between the modes of I and II types. However, one should note that the frequency dependence 
of the critical wave number q c in the one-mode approximation does for I type rolls not show the saturation like in the 
numerical calculations. 

For uoTd ^> 1/ '-\/Pp, when l/(e v ed) H> 1, one can expand Eqs. (|l7| ) in l/(e v Cd)- For the threshold of II type rolls one 
gets 



a» = 



2 (effi + l) m lp 



m 2 h h 



1/2 



— • (19) 



Since <C 1 for cord S> 1 one can see that a c is inversely proportional to wt^. 

By increasing wTd the critical amplitude a c (ujTd) of mode I appears to tend to a constant value already at rather low 
Lord, where e" 1 <C 1. Such behavior was found in [pij by an one- mode analysis. One can show that this feature is not 
a artefact of the one-time-mode approximation. The same result can be obtained by inspection of the high-frequency 
limit in (see Appendix ||) . 

Summarizing, we conclude from our analysis that increasing of the oscillatory flow frequency (for given cell thickness 
d) leads to switching from I to II symmetry type roll instability. It has been shown that this transition is caused by the 
inertia of the nematic fluid described by the parameter eT 1 . For MBBA material parameters one has the transition at 
LJTd ~ 4.1 • 10 4 . Consequently, in the PG experiment flG] with the same liquid crystal and d = 240 /zm the transition 
frequency should occur of about 7 Hz (maximal frequency in this experiment was ~ 2 Hz). Recent experiments [ jl2"| , 
[nil show that this transition could be observed. 



C. Comparison with experiments 

Experiments of Mullin and Peacock (MP) |L2| performed with the nematic E7 show the existence of a roll instability 
at least up to 160 Hz with d between 20 and 50/^m. MP report, that the observation of probe particle have not shown 
any macroscopic flow, which could mean that the rolls are of II, i.e., the oscillating type. We have found in the 
literature only the values of the elastic constants for E7 p3|| , so that we are not able to perform a direct comparison 
of our theoretical results with the MP experimental data. One can consider the influence of material parameters 
variation on the instability type. The results of this consideration are presented in table [n| A good test for the 
quality of oscillatory Couette flow realized in experiments is scaling law, which says, that A/d = F(ujt^), where F is 
some function to be determined from nematodynamic equations. The experimental data of MP |12j satisfy this test 
only in a rather narrow region of low frequencies, namely, up to about 40 Hz. At higher frequencies the scaling law 
becomes broken, which could be caused by failure of surface anchoring. 

Oscillatory shear flow was investigated experimentally also by Anikeev and Kapustina | |l3fl . The authors used an 
eutectic mixture of MBBA and EBBA, which is known to have almost the same material parameters as MBBA. A 
geometry where the director is anchored in the plane of the layer at an angle ip to the shear direction was investigated. 
For the case ip = 90°, which we are interested in, the frequency dependence of a c and q c are presented. The 
experimental points taken from [fl3f are given in Fig.|l| together with theoretical curves. The agreement between 
experimental and theoretical data is good. 



D. Influence of electric field on the instability type and symmetry 

Applying an electric field across the nematic layer can lead to the switching in particular between the different types 
of temporal symmetry of rolls induced by oscillatory shear flow. In the case of MBBA a change of the roll instability 
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of I type to II type by increasing the electric field strength was observed [[[(J . We performed the stability analysis 
of the linearized equations (H)-(||) in a wide range of oscillatory flow frequency and electric field strength for MBBA, 
which has a negative anisotropy of the dielectric permeativity. We have found that the electric field can change not 
only the temporal symmetry of instability, but also the spatial one. The phase diagram (Fig. ||) presents the area, 
where the rolls with symmetry type III appear at threshold. According to their spatial symmetry these rolls represent 
a double-layered structure, so that the critical wave number in this area is about two times larger than otherwise. For 
low values of ojTd an applied electric field leads to an increase of the area of homogeneous instability (HI), so that 
the electric field suppress the roll instability. The profiles of the velocity and the director in the various modes are 
plotted in Fig. ||. 

We have also calculated the phase diagram for a nematic with positive dielectric anisotropy (Fig.^j). The material 
parameters were taken from MBBA except for setting the dielectric anisotropy to be positive. In this case the phase 
diagram (Figji]) presents an extensive area of rolls of II symmetry type. At low uiTd there is a switching between 
homogeneous instabilities of I and II types by varying the electric field. Note, that under these conditions Freedericksz 
transition takes place at voltage Uf (~ 3.74 V). As a result, when U approaches Uf, the critical shear amplitude a c 
goes to zero. 

The one-mode analysis in this case gives unsatisfactory agreement with the full numerical calculations, although it 
keeps all the important features (switching between modes etc.). 



IV. DISCUSSION AND CONCLUSION 



We obtained the complete phase diagram of the orientational instabilities, both homogeneous and spatially periodic, 
in nematic liquid crystal subjected to oscillatory Couette flow based on the linear stability analysis of the full set of 
nematodynamic equations. It has been found for the first time that the inertia of nematic fluid, usually neglected in 
the theoretical considerations, is responsible for the new type of roll instability and switching between roll instabilities 
with different time-symmetry. The results of approximate analytical analysis show an essential features of the influence 
of inertia of nematic fluid on the instability scenarios and are in a good agreement with full numerical simulations 
of underlying equations. For the nematic liquid crystal MBBA increasing of flow frequency leads to the transition 
between rolls of I and II symmetry types. 

Both these modes can be easily distinguished in the experiments. Although there are no reports concerning the ob- 
servations of this transition, the recent experiments on oscillatory flow in nematics show that the technical possibilities 
are rich enough to observe such effect. 

We have also investigated the influence of additional electric field applied across the nematic layer on the flow 
instabilities for the liquid crystals with negative and positive dielectric anisotropy. In case of MBBA it was found that 
the electric field causes transition to the double-layered roll structure, which is intermediate regime between rolls of 
I and II types observed in |10|. Another interesting point in the experimental investigations could be the transition 
between rolls of the same symmetry but quite different wave numbers at increasing of flow frequency (i?7IIi and -R/II2 
in Fig.2). 
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Coefficients for are 



APPENDIX A: COEFFICIENTS 



mi = (ir 2 k 2 + q 2 k 3 )(q 2 r]' 2 + ir 2 )/ri, 

m 2 = (l-ArV^-^^+T^/n, 

m 3 = (1 - A) 7n, 

h = r 2 [7T 2 + k 3 (q 2 -sgn(e a )7r 2 E 2 )]/r 3 , ( A1 ) 

h = (l-AJ-^Ars 



a„-a, 
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with 



n = « 2 (^ + (1-A)-X)+t 3 , 

r 2 = g 2 [7r 2 (ai+ryi+^)+g 2 ^]+7T 4 ryi, 

r 3 = r 2 + (l-A)- 1 [ g X-9 2 ^ 2 (M + «3) + ^ 4 Aa! i ]. 



APPENDIX B: HIGH-FREQUENCY LIMIT 



To investigate the behavior of the critical amplitude in the high-frequency range we neglect e^ 1 in (|1J) and (15) and 
expand the functions / and g in small parameter td = l/(u)Td)- Formally this analysis is valid for 1 <C ujt^ <C 1/P P 
(for MBBA material parameters and d = 20 jam this corresponds to 0.02 Hz <C / <C 10 6 Hz). 

Setting e" 1 = in (|l^)-(|l5|) one has in matrix form 



dt 



Y = A(t)Y, 



(Bl) 



where 



A = 



— edWi —1712a cost 
— ha cost —£dh 



We seek periodic solution of (Bl) of the form 



Y = Y + e^Yx 



At order in e° d one has 



Y = A (i)Y 



(B2) 



where Aq = A 



Since Aq(£) commutes with Ao(i') (unlike A) the solution of (B2) has a form 



Y = eJ 



J^drA (r) Y a Y o = ( f° 



(B3) 



where f° and g° are initial values (fluctuations) of the functions /o and go, respectively. It can be easily seen, that 
Yo is periodic with time averages (fo)t ~ f°, (go)t ~ <7°, so that the state with (f° = 0, g° ^ 0) corresponds to the 
mode of I symmetry type, the state (/° ^ 0, g° — 0) corresponds to mode of II symmetry type. At first order in td 
one has the system 



^I=A (t)Y 1 + b, with b=(-^/° 
dt V — *iffo 



The solution of (B4) is given by 

Yi = Y + J> Ao(t) f dre- & d S A °&b(r). 
Jo 

Substituting expressions for b into ( [B5| ) and using the relation pp| | 

oo 

Io(z) + 2 ^(-) fc / 2fe+ i(z) sin(2fc + 1)9 + 

OO 

2Y / (-) k l2k(z)cos2k8, 



(B4) 



(B5) 



zsinO 



k=Q 



fc=0 



where are modified Bessel functions of the first kind, one has the periodicity condition of Yi 
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io(2\An 2 /2a) 



_ h+mi 
li —mi 

_ h+mi 



for I type, 
for II type. 



(B6) 
(B7) 



Eqs. (B6) and (B7) represent the connection between a and q in implicit form, i.e. the neutral curve ao(q). Minimizing 
ao(q) gives the critical values of the flow amplitudes a c and the wave number q c . Equation (B7) has no real roots for 
q > 0.5, so that when the fluid inertia term is neglected only a roll instability of I symmetry type or homogeneous 
instability of II symmetry type can appear. Eq. ( |B6| ) gives a critical ampli tud e a c ~ 2.2 and wave number q c ~ 12 
and for , that is in a good agreement with numerical results (Fig.|]). From (B7) it follows that the smallest threshold 
a c = 10.9 corresponds to q c = 0. 

The above analysis shows the role of the elasticity described by terms proportional to in the nematodynamic 
equations. The nature of these terms is singular, i.e. when neglecting these terms the critical flow amplitude vanishes 
(a c = 0), although for infinitesimal there exists a finite value of a c . Another point is that at high enough frequencies 
a c does not depend on the flow frequency, which is in contrast to the case of director alignment in the flow plane [0, 
where the critical amplitude of roll formation under oscillatory shear flow decreases with increasing flow frequency. 



APPENDIX C: MATERIAL PARAMETERS 

The following values of the MBBA material parameters at 25 °C pBLEfl: 
K n = 6.66, K 22 = 4.2, K 33 = 8.61 (in units 10~ 12 N); 

ai = -18.1, a 2 = -110.4, a 3 = -1.1, 04 = 82.6, a 5 = 77.9, a 6 = -33.6 (in units 10~ 3 N-s/m 2 ). 
e a = -0.53. 

Freedericksz transition voltage is defined by Uf = K\ K \ xl \ ■ 
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v x 
v z 
n x 



even 
even 
even 
even 



t 

"odd" 

even 
odd 
even 



II 



even 
even 
even 
even 



even 
odd 
even 
odd 



III 



odd 
odd 
odd 
odd 



t 

"odd" 

even 
odd 
even 



IV 



odd 
odd 
odd 
odd 



t 

even 
odd 
even 
odd 



TABLE I. Spatial (2) and temporal (t) symmetries of solutions of linearized equations (|3|)-(|6j). 





Standard values 


Transition values 


New a c 


New q c 


ct\ 1 a-i 


0.164 


1.01 


2.38 


6.9 


ai/ot2 


-0.748 


-5.92 


7.72 


10.1 




-0.706 


-1.83 


6.64 


1.4 


a 6 /a 2 


0.304 


-0.85 


6.60 


1.4 




1.586 


5.95 


2.93 


7.3 


K33 1 K22 


2.05 


0.15 


1.66 


16.6 



TABLE II. In the second column the normalized standard values for the material parameters of MBBA at 25°C are given. 
With those values and at Lord = 825 (typical value for experiment [^2j) the roll instability of I type (overturning rolls) takes 
place with critical amplitude a c = 2.6 and critical wave number q c = 7. When any one of the parameters is varied alone a 
transition to II type rolls (oscillating ones) occurs at the value in the third column. In the fourth and fifth columns the new 
values of the threshold amplitude and the critical wave number are given. 
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FIG. 1. Critical amplitude a c vs. uit^ for MBBA. Solid lines are numerical calculations, dashed lines are approximations 
([t7|), points are experimental data from p3| . RI denotes roll instability, HI is homogeneous instability, Roman numbers denote 
the group of symmetry from table |l| 
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FIG. 2. Phase diagram of instabilities induced by oscillatory shear flow and affected by electric field a nematic with negative 
dielectric anisotropy. RI denotes roll instability, HI is homogeneous instability, Roman numbers denote the group of symmetry 
from table fl. 
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FIG. 3. Director (left) and velocity (right) profiles: 

a) RI I, time averaged, 

b) fl/IIi, t = 0, 

c) i?/II 2 , t = 0, 

d) RI III, time averaged. 
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FIG. 4. Phase diagram of instabilities induced by oscillatory shear flow and affected by electric field in NLC with positive 
dielectric anisotropy. Uf is Freedericksz transition voltage. RI denotes roll instability, HI is homogeneous instability, Roman 
numbers denote the group of symmetry from table F[ 
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